Noise and Controllability: suppression of controllability in large quantum systems 
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A closed quantum system is defined as completely controllable if an arbitrary unitary transfor- 
mation can be executed using the available controls. In practice, control fields are a source of 
unavoidable noise. Can one design control fields such that the effect of noise is negligible on the 
time-scale of the transformation? Complete controllability in practice requires that the effect of 
noise can be suppressed for an arbitrary transformation. The present study considers a paradigm of 
control, where the Lie-algebraic structure of the control Hamiltonian is fixed, while the size of the 
system increases, determined by the dimension of the Hilbert space representation of the algebra. 
We show that for large quantum systems, generic noise in the controls dominates for a typical class 
of target transformations i.e., complete controllability is destroyed by the noise. 
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Coherent control was constructed to steer a quantum 
system from an initial state to a target state via an ex- 
ternal field [l|, • The idea is to control the interference 
pathway governing the dynamics. For pure initial and 
final (target) states the method can be termed state-to- 
state coherent control. A generalization is steering si- 
multaneously a set of initial pure states to a set of final 
states, i.e. controlling a unitary transformation. Such an 
application sets the foundation for a quantum gate oper- 
ation [^-[5|. Three basic questions address feasibility of 
coherent control. The first, for a preset initial and target 
state does a control filed exist? This is the problem of 
controllability. The second, how to construct a field that 
leads to the target? This is the problem of synthesis. The 
third, how to optimize the field that carries out this task? 
This is the problem of Optimal Control Theory [^-Q. 
Experimentally there has been a remarkable success in 
constructing devices designed to generate arbitrary con- 
trol fields |9l-[ll|. Nevertheless, in practice controllability 
is hard to achieve even for small quantum systems |12| - 
Applications toward quantum information process- 
ing require upscaling of the control procedures to large 
quantum systems. The issue of controllability of a closed 
quantum system has been addressed by Tarn and Clark 
[15{ . Their theorem states that for a finite dimensional 
closed quantum system, if the control operators and the 
unperturbed Hamiltonian generate the Lie-algebra of all 
Hermitian operators, the system is completely control- 
lable, i.e., an arbitrary unitary transformation of the sys- 
tem can be realized by an appropriate application of the 
controls [l^. Complete controllability implies state-to- 
state controllability. 

In practice the controlled systems are open and the 
system-bath coupling introduces noise into the system 
dynamics. A number of techniques have been designed 
to combat the environmental noise and effective and in- 
genious algorithms have been invented and explored [17- 
UV{ . The present study focuses on the effect of the noise 
originating in the control field. The magnitude of this 



noise depends on the properties of the control field. This 
dependence raises a fundamental question: is it always 
possible, for a given target, to design a field, such that 
the effect of the associated noise can be neglected? This 
problem has been extensively investigated in the context 
of fault-tolerant quantum computation [2^ , where vari- 
ous schemes have been designed to fight the noise in the 
gates. In quantum computation the number of gates in- 
creases with the size of the system. In many fields, for 
example in NMR [235 o'" control of molecular systems 
1] a different control paradigm is standard. There the 
control operators are fixed while the size of the system 
may vary. This is the control paradigm considered in the 
present study. The presentation is restricted to quan- 
tum system with finite Hilbert space dimension. It is 
assumed that the control operators are elements of the 
spectrum-generating algebra of the quantum system [23] . 
For a finite dimensional system it is sufficient to consider 
a compact semisimple algebra [1^ . The size of the sys- 
tem is determined by the dimension of the Hilbert-space 
representation of the algebra. It can be characterized by 
a parameter N , which is a scaling factor of the highest 
weight of the representation [2^ . The physical interpre- 
tation of N depends on the system. It can be, for exam- 
ple, a number of particles in the system or the number 
of energy levels. For this paradigm it is shown that com- 
plete controllability of a large quantum system, i.e., for 
^ 1, even in the weaker sense of state to state control, 
does not survive in the presence of generic noise in the 
control fields. As a result the control is not scalable with 
the size of the system. 

Let the Hamiltonian of the controlled system be 



H = H, 







^K(i) + a(t)]x, 



(1) 



where Ufe(t) are control fields and the noise S,k{t) is 
delta-correlated Gaussian noise with {S,k{t)^i{t')) = 
2ri\uk{t)\SkiS{t — t'). The dimensionless number rj > 
measures the relative strength of noise in the control 
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fields. The operators are elements of the spectrum 
generating algebra associated with Hq and are termed 
controls in what follows. The equation of motion for the 
density operator of the system is given by [13] 
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In the absence of noise the system ([2]) is completely 
controllable for a generic Hq. Due to the noise the purity 
V = Tr|/5^} of an initially pure state p = \ip) {ip\ will 
decrease. The effect of noise can be neglected only if 
the purity loss during the target transformation is small, 
i.e., AV <C 1. Otherwise, the state-to-state and, hence, 
complete controllability is lost. The instantaneous rate of 
purity loss for a pure state p — lip) {ip\ evolving according 
to Eq.© is given by 



dt 



4, E 



where 
in the state ip- 



(3) 

is the variance of the control operator Xfc 



Xt 



(4) 



It should be noted that the rate of purity loss increases 
with the variance of X^ . A generic state of the system is 
characterized by Aj^^ [tp] ~ N^. 

The logic of the presentation is as follows. First a 
class of state-to-state target transformations is defined. 
The transformations can be accomplished in the absence 
of noise due to complete controllability. The time dura- 
tion of these transformations can be bounded from below. 
Next the lower bound of the purity loss rate for the evolv- 
ing state in the presence of noise is estimated. For small 
purity loss the two bounds can be combined to obtain the 
lower bound on the purity loss during the transformation. 
This bound depends on the relative strength of noise rj. 
The crux of the class of transformations considered is 
that for any realization of the control fields accomplish- 
ing the transformation the evolving system resides for a 
long period of time in states with large variance ~ N"^ 
with respect to the control operators. The large variance 
will generate ~ ijN'^ rate of purity loss in the presence 
of noise in the control fields. For this class of transfor- 
mations it is found that the lower bound on the time 
of transformation scales as ~ N~^, therefore, r] must 
be 0{N~^) in order that the purity loss be negligible. 
Since in practice the relative magnitude of the noise can- 
not be made arbitrarily small, it follows that the loss of 
purity cannot be neglected for large systems. As a re- 
sult, a transformation from an arbitrary pure state to 
another pure state cannot be accomplished. Large sys- 
tems are not state-to-state controllable and are therefore 
not completely controllable. 



Let |n) denote the eigenstates of Hq. We 
consider the transformation from the initial state 
\tpi) = I"-) to the final state IV'/) — 

J2n'^f,ne"^^'" \n), such that 



1 > IIArll > £ > 0, 



(5) 



where ||Ar|| is the Euclidean norm of Ar = r/ — r^, 
i"i iri,i,ri^2, ■■■) and r/ = (r/,i, r/,2, ■■■)■ The choice 
of norm excludes changes to states that can be reached 
by free propagation generated by Hq- As a consequence, 
the free evolution does not contribute to the bound on 
the time of the transformation ipi ^ ipf. Under assump- 
tion that the noise is small the bound can be estimated in 
the zero order in the noise strength. For estimation of the 
bound an auxiliary operator A is defined such that: (i) 
it commutes with Hq; (ii) its expectation value changes 
during the transformation. Since A commutes with Hq 
the change of its expectation value during the transfor- 
mation is due to the operation of the control fields. We 
define A = J2n !"■) ("■!' where s„ — sign{Ar„}. The 
change of the expectation value of the operator A during 
the transformation ipi ^ ipf is given by 
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where we have used the fact that the vector of ampli- 
tudes r is nonnegative, and, therefore, |Ar„| > Vi^n only 
if Ar„ > 0. Using inequality ([5]) we obtain 



A 



A >e' 



(7) 



which gives the minimal change of the expectation value 
of the operator A during the transformation ipi — ipf. 
On the other hand, the change of the expectation value 
of A can be estimated from the Heisenberg equations: 



dt 



i^Uk{t) 



where we have used the fact that 
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Ho, A 



(8) 



0. Let the 



time of the transformation be T. Then, 
'A^ - (A) ^ ^^(A) dt 

T 



< > / dt max 

L 0<t<T 



Xfc, A 



< 



2E r\Mt)\ dt\Ak\, 
u Jo 



(9) 



where A^ ~ A?" stands for the eigenvalue of X^, maximal 
by the absolute value. In the derivation we have used the 



3 



fact that the eigenvalues of A are ±1. Defining the av- 
erage control amplitude Uk = ^ \uk{t)\ dt, and using 
Eqs.dll) and we arrive at the inequality 



T> e'{2Ek^k\Ak\)-'-eH2Nj:, 



Uk 



(10) 



which bounds the time of the transformation for given 
Uk- This bound is similar to bounds obtained for the 
transformation to an orthogonal state in Refs.[29l|. 

Without loss of generality we assume that the time 
T of the transformation 7/;^ — >■ -0/ is the first passage 
time when the evolving state \ip{t)) = J2n ''r! (i)e*'^"^*^ \n) 
"crosses the border" ||r(T) — ri|| ~ e, i.e., such that 
for t < T we have \\r{t) — r^H < e. Under assump- 
tion that the purity loss AV during the transforma- 
tion is small, the evolving state can be approximated 
by p{t) = -f- w = iV'(t)) (V^i)!- Taking the 
leading contribution of p'"'^' into account, we estimate the 
lower bound on the purity loss from Eq.Q: 

AV > 4TrjJ2uk ^nnn^{A^jm] 



^(i) [Xfe,[Xfe,p(i)(t)]] m)}- (11 



We further assume that during the transformation 
^Xfc bPi^)] ~ (^fe)^ ^ III this case we can neglect 
the pf^^-dependent term in the inequality (|lip . Using 
the inequality ([TU| . we obtain 

AV > ^. (12) 



To estimate mino<(<T j^Aj^^ [V'(i)]| we find the lower 



bound on the variance of Xfc in the states \ip) — 
^„r„e*'^"|n) such that ||r — < e. The variance 
^Xfc t^''] ^ function of the amplitudes r = (ri,r2,...) 
and the phases (t>i, (f)2, ■■■■ The free evolution can change 
the phases at no cost in purity. Therefore, the minimal 
variance attainable for given amplitudes is sought: 



Aj^ (r) = min 



(13) 



We assume, that Aj^^ (r) is a smooth function of r for 
I! r — Till < ||Ar||, i.e., for sufficiently small Ar and 
\Sr\ < ||Ar|| we can expand Aj^^ (r^ + Sr) « Aj^^ (r^) + 
VAj^ (r^) • Sr. We note that 



S^x. (r) ^ VA.^ (r) • Sv 



< 



VAx.. (r) 



\Sr\\. (14) 



Let the minimum in Eq. (fT3)) obtain at 05!(r), 02(r), 
Let us denote the associated state as ■i/'*- Then A-j^ (r) 
Ax, [^1 and 



VAx^_ (r) = WA. [r] 



(15) 



It should be noted that Ay^ [ij;*] depends on r„ both 
through the amplitudes of t^* and through the phases 
0*(r), which are also functions of r. Nonetheless, since 
0*(r) are defined as giving the minimum of Aj^^ [ip], 
derivatives of Aj^^ [tfj*] with respect to 4>n{^) vanish and 
0*(r) may be considered as r- independent for the oper- 
ator V in the rhs of Eg. dTS]) . 

Using the definition Eq.Q we obtain 



VAx fr) = yU 



2U* Xk r)y{r ^ 



(16) 



Using the explicit form of ip* and the fact that V act 
only on the state's amplitudes we can show that the 
Euclidean norm of the rhs of Ea. ([TB)) is bounded by 

3%/2(Afc)^ Then, Eqs. dH]) and (HH) imply 5A^^ (r) 

3\/2(Afc)^ ||(5r||. It follows that for \^) = X]„r„e^'^" 
with llr - r;|| < e < 1 



< 



A^jV']>AxJr,)-3V2(Afc)^e (17) 
From inequalities (IT2|) and p7)) we obtain 

/minJAx,(r,)-3^/2(A0'e|\ 

> V ^ ■ TTTn (18) 

I max/{|A/|} I 

The variance Aj^^ [ip] scales as iV^ in a generic state of the 

system . The scaling of Ay^^ (r^) is a more subtle ques- 
tion, since it is the outcome of the minimization with 
respect to the phases in the eigenstates basis. For our 
purposes it is sufficient to show that A-j^^ (r^) ~ N'^ for 

some Let's consider a generic eigenstate \(f>) of Hg. 
The variance Ay^ [(/>] scales as N'^ for all k. Moreover, 
the variance is independent of phases. Therefore, tak- 
ing j-f/^i) — \(f>) we shall have Aj^^ (r^) ~ iV^, and, for 
sufficiently small e, inequality (jlSp will imply 

AV > 2e^r]N. (19) 

It is important to note that inequality ([T5)) holds for small 
e, which can still be of the order of unity with respect 
to N. As a consequence, for such target transformations 
the rhs of inequality (|19p will scale as 77 A^. In order that 
the purity loss be negligible, the relative noise strength 77 
must scale as 0{N~^), which is unrealistic for large N. 

As an example, let's consider a system of cold atoms in 
a double well trap [lol ■ The Hamiltonian of the system 
can be put in the form 



Hn 



-LoJx + S3z + 



A^' 



(20) 



where Jfc are elements of the su{2) Lie-algebra. The oper- 
ator Jz corresponds to the population difference between 
the wells and 3y to the population flow between the wells. 
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The on-site interaction strength is given by U, w is the 
hopping rate of the atoms and 6 determines the tilt of the 
potential wells. N is the number of atoms, corresponding 
to the iV + 1 dimensional irreducible representation of the 
su{2) algebra, with the spin number j = iV/2. Control 
can be attained via the control operators Ja; and J2, i.e., 
by modulating the hopping rate and the tilt of the wells 
[31l] . For a typical system of TV > 10^ atoms the nec- 
essary condition rjN ^ 1 cannot be satisfied in practice. 
Therefore, the purity loss in the system during a target 
transformation of the type considered in the present work 
will be of the order of unity and the effect of noise cannot 
be neglected. 

We conclude that the statc-to-state controllability of 
large quantum systems is destroyed by the noise on the 
control. The strategy of suppressing the influence of the 
noise by a faster control will fail due to the required in- 
crease in field amplitude, inevitably accompanied by the 
increase in noise. The purity loss of the evolving state 
will be faster due to the increase in noise on the con- 
trols. The order-of-unity decrease of the purity implies 



that the relative error in the expectation values of some 
of the system operators in the target state is of the or- 
der of unity. This large relative error is characteristic 
for transformations between states, where variance of the 
control operators is ^ iV^. It is expected that the rela- 
tive error is small for target transformations between the 
"classical" generalized coherent states [32l| with respect 
to the spectrum generating algebra, where the maximal 
variance of the control operators is ~ TV. An exarnple is 
transformations between the spin-coherent states [3J| of 
a quantum spin, corresponding to BEC states of atoms 
in the double-well trap. Finally, it is interesting to note 
that the scalability of the control paradigm considered in 
the present study can be attacked from a different angle 
of resource management in quantum computing fsi*] . Re- 
markably, the notion of what constitutes a large quantum 
system is similar in both studies. 

Acknowledgements We are grateful to Kalvi institute 
University of Santa Barbara for hospitality. Work sup- 
ported by the Israeli Science foundation. 



S. Rice, Science 258, 412 (1992). [18 
M. Shapiro and P. Bmmer, Rep. Prog. Phys. 66, 859 
(2003). [19 
C. M. Tesch, L. Kurtz and R. de Vivie-Riedle, Chem. 
Phys. Lett. 343, 633 (2001). [20 
J. P. Palao and R. Kosloff, Phys. Rev. Lett. 89, 188301 

(2002) . [21 
H. Rabitz, M. Hsieh, and C. Rosenthal, Phys. Rev. A 72, 
052337 (2005). [22 
S. H. Shi, A. Woody, and H. Rabitz, J. Chem. Phys. 88, 
6870 (1988). [23 
R. Kosloff, S. A. Rice, P. Gaspard, S. Tersigni and D. 
Tannor, Chem. Phys. 139, 201 (1989). [24 
J. P. Palao and R. Kosloff, Phys. Rev. A 68, 062308 

(2003) . 

T. Weinacht, J. Ahn, and P. Bucksbaum, Nature 397, [25 
233 (1999). 

M. Aeschlimann, M. Bauer, D. Bayer, T. Brixner, F. J. [26 
Garcia de Abajo, W. Pfeiffer, M. Rohmer, C. Spindler, [27 
and F. Steeb, Nature 446, 301 (2007). 

R. Selle, P. Nuernberger, F. Langhojer, F. Dimler, [28 
S. Fechner, G. Gerber, and T. Brixner, Opt. Lett. 33, [29 
803 (2008). 

A. Walther, B. Julsgaard, L. Rippe, Y. Ying, S. KroU, 
R. Fisher, and S. Glaser, Physica Scripta T137 (2009). [30 
A. Sporl, T. Schulte-Herbruggen, S. J. Glaser, 
V. Bergholm, M. J. Storcz, J. Ferber, and F. K. Wil- [31 
helm, Phys. Rev. A 75, 012302 (2007). 

S. Chaudhury, S. Merkel, T. Herr, A. Silberfarb, I. H. [32 
Deutsch and P. S. Jessen, Phys. Rev. Lett. 99, 163002 
(2007). [33 
J. Clark and T. Tarn, J. Math. Phys. 24, 2608 (1983). 
V. Ramakrishna and H. Rabitz, J. Math. Phys. 54, 1715 [34 
(1996). 

[17] L. Viola, E. KniU and S. Lloyd, Phys. Rev. Lett. 82, 2417 
(1999). 



K. Khodjasteh and L. Viola, Phys. Rev. Lett. 102, 
080501 (2009). 

K. Khodjasteh and D. A. Lidar, Phys. Rev. Lett. 95, 
180501 (2005). 

K. Khodjasteh and D. A. Lidar, Phys. Rev. A 75, 062310 
(2007). 

G. Gordon, N. Erez and G. Kurizki, J. Phys B: At. Mol. 
Phys. 40, S75 (2007). 

J. Preskill, Proc. R. Soc. Lond. A 454, 385 (1998), 

E. KniU, Nature. 434, 39 (2005). 

I. I. Maximov J. Salomon G. Turinci and N. C. Nielsen, 
J. Chem. Phys. 132, 084107 (2010). 

A. Bohm, Y. Neeman, and A. O. Barut, Dynamical 
groups and spectrum generating algebras (World Scien- 
tific, Singapur, 1988). 

R. Gilmore, Lie groups, Lie Algebras and Some of Their 

Applications (John Wiley and Sons, New York, 1974). 

R. Gilmore, J. Math. Phys. 20, 891 (1979). 

V. Gorini and A. Kossakowski, J. Math. Phys. 17, 1298 

(1976). 

S. Boixo, L. Viola, and G. Ortiz, EPL 79, 40003 (2007). 
N. Margolus and L. B. Levitin, Physica (Amsterdam) 
120D, 188 (1998), L. Levitin, T. ToffoU, and Z. Walton, 
Int. J. Theor. Phys. 44, 965 (2005). 

J. R. Anglin and A. Vardi, Phys. Rev. A 64, 013605 
(2001). 

S. Levy, E. Lahoud, and I. Shomroni, J. Steinhauer, Na- 
ture 449, 579 (2007). 

W. Zhang, D. H. Feng, and R. Gilmore, Rev. Mod. Phys. 
62, 867 (1990). 

F. Arecchi, E. Courtens, R. Gilmore, and H. Thomas, 
Phys. Rev. A 6, 2211 (1972). 

R. Blume-Kohout, C. Caves, and I. Deutsch, Found. 
Phys. 32, 1641 (2002). 



